Two-photon quantum interference from separate nitrogen vacancy centers in diamond 
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We report on the observation of quantum interference of the emission from two separate nitro- 
gen vacancy (NV) centers in diamond. Taking advantage of optically induced spin polarization in 
combination with polarization filtering, we isolate a single transition within the zero-phonon line 
of the nonresonantly excited NV centers. The time-resolved two-photon interference contrast of 
this filtered emission reaches 66%. Furthermore, we observe quantum interference from dissimilar 
NV centers tuned into resonance through the dc Stark effect. These results pave the way towards 
measurement-based entanglement between remote NV centers and the realization of quantum net- 
works with solid-state spins. 



The nitrogen vacancy (NV) center in diamond is an at- 
tractive candidate for the construction of quantum net- 
works [T]-[3]. It is a highly stable single photon emitter 
[4], whose optical transitions are linked [5 to a long- 
lived electronic spin with excellent coherence properties 
[6]. The electronic spin can be coherently manipulated 
with high fidelity with pulsed microwave fields ^ . By 
extending these techniques to proximal nuclear spins, a 
controllable multiqubit quantum register can be realized 

[Ml]. 

Optical schemes based on photon interference provide 
a way to establish quantum interactions over a distance. 
When two photons that are indistinguishable simultane- 
ously impinge on a beamsplitter they coalesce into the 
same output port [12 . If these two photons are each 
entangled with the spin of the emitter, their interfer- 
ence can be exploited to generate entanglement between 
two distant emitters p^SHTS] . This type of measurement- 
based entanglement has recently been achieved between 
two remote single ions [16] . While spin-photon entangle- 
ment has been realized using NV centers [5 , two-photon 
quantum interference (TPQI) with NV centers has so far 
not been observed. 

Here we capitalize on recent advances in the un- 
derstanding of the NV center's excited states [TTHII] 
and a drastic improvement of the photon collection effi- 
ciency [ITII22], and demonstrate TPQI from two separate 
NV centers. Furthermore, we show that even dissimilar 
NV centers can exhibit TPQI when their emission lines 
are brought into resonance using dc Stark tuning [211 [23] . 
This latter technique greatly facilitates a scaling to larger 
networks. 

We perform our experiments on naturally occurring 
NV centers in high purity type Ila chemical-vapor depo- 
sition grown diamond [24] with a (111) crystal orientation 
obtained by cleaving a (100) substrate. Because a TPQI 



measurement involves coincidence detection of photons 
emitted from two centers, high collection efficiency is es- 
sential. To this end, we deterministically fabricate solid 
immersion lenses ^ around preselected centers by 
focused ion beam milling. Figure 1(a) shows an optical 
microscope image of the device being used; the inset is an 
scanning electron microscope image of a similar device. 
In the confocal scans [Fig. 1(b)], the NV centers appear 
as bright spots inside the lenses. To enable spin manipu- 
lation and Stark tuning of the optical transition energies, 
we lithographically define a gold microwave stripline and 
gates around the solid immersion lenses. A dual path 
confocal setup allows us to individually address two NV 
centers within the same diamond [Fig. 1(c)]. The sample 
is mounted inside a flow cryostat and experiments are 
performed at a temperature of 9K. 

The optical emission spectrum of the NV center 
[Fig. 2(a)] consists of both direct transitions between the 
ground and the excited state [the so-called zero phonon 
line, ZPL, that contributes 4% to the total emission] 
and transitions that additionally involve the emission of 
phonons [phonon side band (PSB)]. At low temperatures, 
the ZPL emission spectrum exhibits several narrow lines 
which, for low-strain centers, correspond to spin selective 
transitions between the ground and excited state [T7] . 
Observation of TPQI requires indistinguishable photons 
which we produce by isolating a single transition within 
the ZPL. Appropriate bandpass filters remove the inco- 
herent fraction of the emission (PSB). Nonresonant exci- 
tation, as used in this experiment, polarizes the NV elec- 
tronic spin into the = state, hence only transitions 
between the and Ey excited states (m^ = 0) and the 
= ground state level occur. The dipoles associated 
with these two transitions are orthogonal to each other 
and to the N-V axis; by working with NV centers ori- 
ented along the (111) direction we ensure that collected 
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FIG. 1. (a) Optical microscope image of the sample. The 
microwave stripline (MW) is used for spin manipulation and 
the gates (Gl and G2) are used for Stark tuning of the opti- 
cal transition energies. Inset: Scanning electron microscope 
image of a similar device, (b) Confocal microscope image 
showing NV A and NV B (logarithmic color scale), (c) Exper- 
imental setup. The two NV centers inside the same diamond 
are simultaneously excited with either a resonant (637 nm) or 
an nonresonant (532 nm) laser. Two separate paths allow de- 
tection of photons emitted into the zero-phonon line (ZPL) 
or the phonon side band (PSB). A variable retardance wave- 
plate (WP) aligns the polarization of the desired transition 
to the polarizing beamsplitter (PBS) while compensating for 
ellipticity introduced by other optical components. 

ZPL photons remain linearly polarized, with orthogonal 
polarizations for Ex and Ey |20] . Consequently, for 
a center that is spin polarized into mg = 0, we can iso- 
late the Ex or Ey emission line by placing a polarizer in 
the detection path. Furthermore, we maximize the signal 
by using the polarization selectivity to excite one dipole 
predominantly. 

To characterize the ZPL fine structure, we use pho- 
toluminescence excitation (PLE) to probe the ZPL ab- 
sorption spectrum and a scanning Fabry- Perot (FP) cav- 
ity to analyze the ZPL emission spectrum. PLE spectra 
are obtained by sweeping the frequency of the resonant 
excitation laser [637 nm. Fig. 1(c)] across the ZPL tran- 
sitions while recording the red-shifted emission into the 
PSB. During PLE scans, we ensure the correct charge 
state and visibility of the = transitions by simul- 
taneous weak green excitation [26^. Figure 2(b) shows 
the PLE spectra of the two selected emitters. While the 
lower-energy Ey transitions are 2.9 GHz apart from each 
other, the higher-energy Ex transitions partly overlap. 
This near overlap is also observed by using a scanning FP 
cavity to monitor the ZPL emission under nonresonant 
(532 nm) excitation [Fig. 2(c)]. By setting the filters to 




FIG. 2. Spectral characterization of NV A and NV B. (a) 
Level structure. The rris = ground state is connected via 
spin selective transitions to the Ex and Ey {rris — 0) excited 
state, TTis = ±1 states are omitted for clarity. Photons that 
are emitted into the phonon side band (PSB) can be spec- 
trally separated from resonant emission. Also shown are the 
resonant (637 nm) and the nonresonant laser (532 nm). (b) 
Resonant excitation of the ZPL and detection of the PSB. 
The frequency of the red laser is swept across the resonances 
while weak green excitation prevents photoionization and op- 
tical pumping (see text). The spectra are averaged over 70 
scans, (c) Scanning Fabry-Perot spectrum under nonreso- 
nant excitation. Only the Ex transition is visible confirming 
sufficient suppression of the Ey transition by polarization fil- 
tering. The hnewidths (450 MHz for NV A and 570 MHz for 
NV B) include drifts of the cavity during the measurement. 
The offset is due to detector dark counts. 



transmit only the linearly polarized light associated with 
the Ex transitions, we observe a single emission line from 
each NV center with nearly equal frequencies. 

Inhomogeneous broadening of the emission Hnewidths 
will limit the observability of quantum interference [27J. 
For individual PLE scans, obtained in the absence of 
green light, we find narrow 38±7MHz (36 ± 5 MHz) 
Hnewidths for NV A (NV B); however, a 532 nm re- 
pump pulse between scans leads to an overall distribu- 
tion of frequencies with an inhomogeneous linewidth of 
263 ±6 MHz for NV A and 483 ±5 MHz for N V B [20l [28] . 
The same spectral diffusion is observed in PLE spectra 
recorded with simultaneous green excitation [Fig. 2(b)] 
and in the ZPL emission spectrum under 532 nm excita- 
tion [Fig. 2(c)]. Although the observed inhomogeneous 
broadening exceeds the radiative linewidth by an order 
of magnitude, two-photon interference effects can still be 
detected. Provided that the emission linewidth does not 
exceed the inverse time resolution of the photon detec- 
tors, simultaneous detection of a photon from each NV 
center erases the which-path frequency information, al- 
lowing quantum interference to be observed [29l [30] . 

For the interference measurement, we employ a green 
(532 nm) pulsed laser (62 ps) with a repetition frequency 
of 10 MHz both to excite the two NV centers and to ini- 
tialize the spin state into rris = 0. A combination of 
spectral filtering (ZPL filter) and polarization rejection 
(as discussed above) isolates the Ex lines of NV A and 
B. This emission is then coupled into the two input ports 
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of a polarization-maintaining fiber beamsplitter, ensur- 
ing excellent spatial mode overlap. We establish tem- 
poral mode overlap by using equidistant excitation and 
collection paths for the two centers. The output ports of 
the beam splitter are connected to two avalanche photo 
diodes (APDs) with sub-ns time resolution; the APD 
signals trigger the start and the stop of a fast count- 
ing module whose jitter is less than 12 ps. By recording 
the coincidence counts as a function of start-stop delay, 
we perform a time-resolved measurement. 

As a calibration experiment, we insert a A/2 wave plate 
into one detection path [Fig. 1(c)] so that the two photons 
entering the beamsplitter have orthogonal polarization. 
This makes the two photons distinguishable and no in- 
terference can be observed. Therefore, the coincidence 
distribution only reveals the temporal overlap of two in- 
dependent photon wave packets [see Fig. 3(a)]. The sit- 
uation changes dramatically when the photons enter the 
beamsplitter with parallel polarization [Fig. 3(b)]. For 
this case, two-photon quantum interference is observed: 
around zero detection time difference (r = 0) the two 
photons mainly leave the beamsplitter into the same out- 
put port, leading to a significant reduction in coincidence 
detections. For larger time differences interference is con- 
cealed because of averaging over many photons with dif- 
ferent frequencies [27|. We thereby observe TPQI as a 
reduction in coincidence detection events within a time 
window given by the inverse of the inhomogeneous emis- 
sion linewidth. 

In principle, full visibility interference can be observed 
even for inhomogeneously broadened and detuned pho- 
ton sources at the cost of a reduced width of the in- 
terference dip [29]. The observed contrast at r = is 
limited primarily by NV emission into undesired spectral 
lines owing to imperfect control over the charge and spin 
state of the center. Specifically, when the center is in 
the neutral charge state (NV^), a portion of its broad- 
band fluorescence lies within the 3 nm bandwidth of our 
ZPL filters; from independently measured optical spec- 
tra, these NV^ photons contribute approximately 10% to 
the collected emission. Furthermore, under 532 nm ex- 
citation a residual ~10% = ±1 spin occupation [11] 
produces circularly polarized emission on several other 
transitions, contributing of our polarization-filtered 
signal. From such a 15% background level we expect at 
most 72% visibility. Furthermore, finite time resolution 
in our detection system will average over sharp tempo- 
ral features, raising the depth of the zero-time-difference 
minimum. 

We observe a quantitative agreement between our data 
and no- free-parameter simulations of the experiment. 
Following Legero et al. [27], we model TPQI of exponen- 
tially decaying photon wavepackets with Gaussian fre- 
quency noise and calculate the expected coincidence de- 
tections using independently measured parameters. At 
T = we observe a contrast of 66±10% which is to our 
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FIG. 3. Two-photon quantum interference, (a) Orthogonal 
polarization. Periodic peaks of the coincidence counts cor- 
respond to the 10 MHz repetition frequency of the excita- 
tion. The coincidence distribution shows typical bunching 
and antibunching features of the two independent NV centers 
under pulsed excitation [37]. (b) Parallel polarization. In- 
terference of indistinguishable photons leads to a significant 
decrease in coincidence events at zero time delay. Simula- 
tions (solid line) using only independently measured parame- 
ters are done according to Legero et al. [27 . The dashed lines 
show 1 standard error uncertainty in the simulation. The 
dark and light grey areas illustrate 1 and 2 standard devia- 
tions expected scatter. Parameters used for the simulations 
are: excited state lifetime = 12.0±0.4ns, detuning between 
the two centers = 130±150MHz, frequency jitter between 
them = 550±80MHz, APD jitter = 410=bl0ps, dark count 
rate = 60±10s-\ NV B count rate = 1470=b50s-\ NV A 
count rate = 2700±50s"\ background = 15±5%. 



knowledge the highest value reported for two separate 
solid-state emitters j3Qti33] . This value can be improved 
by more stringent filtering of the ZPL emission or by in- 
creased control over the spin and charge states. We note 
that for measurement-based entanglement, the visibility 
determines the fidelity while the width of the interfer- 
ence dip sets the success probability of the entanglement 
operation. 

To observe TPQI, we selected NV centers with nearly 
identical Er^. frequencies; in general, however, because of 
their different strain environments, two NV centers are 
unlikely to exhibit the same emission frequencies jlSl |34] . 
Even in high purity type Ila samples, we typically observe 
a spread of tens of GHz between different centers. Nev- 
ertheless, NVs can be tuned into resonance by applying 
electrical fields to induce dc Stark shifts of the transition 
energies ^23|. The tuning range of our devices reaches 
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several GHz and is enhanced to tens of GHz in the pres- 
ence of 532 nm excitation [21] . This enables us to observe 
TPQI in the general case of dissimilar sources. 

In the absence of an external voltage bias, the Ey tran- 
sitions of NV A and NV B lie far apart in energy so 
that their TPQI cannot be measured within the time 
resolution of our detectors. By applying a voltage to 
gate 1 [Fig. 1(a)] while keeping the microwave stripline 
and gate 2 on ground, we can tune the Ey transitions into 
resonance. Because the degree of tuning is strongly af- 
fected by the presence of nonresonant excitation [2lj , we 
calibrate the Stark tuning under the same pulsed 532 nm 
excitation conditions as used for the interference mea- 
surement: the PLE spectra in Fig. 4(a) show the ad- 
ditional fluorescence induced by resonant excitation, re- 
vealing the spectral location of the Ey line. Near -13.6 V, 
we observe overlapping transition energies for the two NV 
centers, establishing the appropriate setting for observa- 
tion of TPQI. In addition to shifting the energy of the NV 
center transitions, applied electric fields can also rotate 
the axes of the E^ and Ey dipoles; we find that signifi- 
cantly different polarization settings are required to filter 
the Ey emission at this gate voltage. 

Using the calibrated voltage and polarization settings, 
we measure the time-resolved interference of photons 
emitted on the Ey transitions from each NV center. Fig- 
ure 4(b) shows a significant decrease in coincidence de- 
tection events at zero time difference. The width of the 
interference signal is smaller than what would be ex- 
pected from earlier measurements of the inhomogeneous 
linewidth of the emitters. Because the laser intensity 
changes the effective Stark shift, a drift in sample posi- 
tion directly translates into a spectral shift. To account 
for the spectral variations, we perform the measurement 
by alternating PLE spectroscopy and minute-duration 
coincidence detection, and post-select the data based on 
their relative detunings. We find our data to agree with a 
simulation based on the measured frequency distribution 
[solid line in Fig. 4(b)]. In contrast, a model without 
TPQI [shaded area in Fig. 4(b)] fails to reproduce the 
observed drop in coincidences around r = 0. Increased 
control over spatial and laser power fiuctuations is ex- 
pected to greatly enhance the interference contrast in the 
presence of gate voltage tuning. 

We have demonstrated two-photon quantum interfer- 
ence with separate NV centers in diamond. The observed 
contrast can be further improved through resonant exci- 
tation, which eliminates photon emission from the incor- 
rect charge and spin states. Moreover, the coincidence 
rates may be enhanced by embedding NV centers into 
optical cavities [35[ [36] . When combined with recently- 
demonstrated spin-photon entanglement |5j, our results 
enable remote entanglement of NV centers, and open the 
door to applications in quantum information processing 
and long distance quantum communication. 

We thank B. Hensen, V. Jacques, M. D. Lukin, Y. 
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FIG. 4. Interference of two dissimilar sources, (a) Ey tran- 
sition energy as a function of applied gate voltage. The two 
centers show opposite tuning behavior and are brought into 
resonance at -13.6V. Data are taken under simultaneous green 
excitation with the same power that is used for the interfer- 
ence measurement, (b) Two-photon quantum interference. 
The overall coincidence distribution is the sum of 255 one- 
minute histograms. Before each interference measurement we 
perform a PLE scan. We select only histograms for which 
the relative detuning of the two centers is between 350 and 
1200 MHz. The simulation (red line) is based on the measured 
frequency distribution and shows oscillations that reflect the 
discrete set of detuning values. The dark and light grey areas 
illustrate the expected signal for two noninterfering sources 
with 1 and 2 standard deviations, respectively. 
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